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Abstract 

We give a survey concerning both very classical and recent results on the electro¬ 
static interpretation of the zeros of some well-known families of polynomials, and 
the interplay between these models and the asymptotic distribution of their zeros 
when the degree of the polynomials tends to infinity. The leading role is played by 
the differential equation satisfied by these polynomials. Some new developments, 
applications and open problems are presented. 
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1 Introduction 


The electrostatic interpretation of the zeros of the classical orthogonal polynomials is 
probably one of the most elegant results in the theory of special functions, linked in 
the first instance to such a distinguished name as Stieltjes (although studied also by 
Bocher, Heine, Van Vleck, and Polya). Although this topic has remained “dormant” 
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for almost a century, a renewed interest has appeared recently, partially motivated 
by the connection of this topic with modern powerful techniques from the theory of 
logarithmic potentials, as well as by an increasing interest on the study of new classes 
of special functions with roots in Physics, Combinatorics, Number Theory, etc. 

This paper is a short and light survey on these topics. Our intention is to provide not a 
comprehensive list of results, but more to convey the “flavor” of techniques and ideas 
behind the electrostatic model and its direct implications in the asymptotic theory. 
With this purpose, in Section 2 we derive the classical electrostatic interpretation of 
the zeros of Jacobi polynomials, found originally by Stieltjes, and we describe some 
of its new and still-in-progress direct generalizations (Section 3). Next we discuss a 
recent progress by Ismail in connecting electrostatics with orthogonality in Section 4. 
Seeking for a model for complex zeros in Section 5 we formulate an alternative model, 
not constrained to the real line, as well as discuss some further generalizations and 
applications of the electrostatic model in Section 6. We conclude the paper with some 
examples of applications of these models to the study of the asymptotic distribution 
of zeros in the semiclassical limit. 


2 Electrostatic model for classical orthogonal polynomials 


Let us begin with a classical and very well-known result, found by Stieltjes in 1885, 
about the electrostatic interpretation of the zeros of Jacobi polynomials The 

dehnitions and basic properties of these polynomials can be found in Chapter IV of the 
classical monograph by Szego [63]; we will mention here only those that will be needed 
further. 


Jacobi polynomials can be given explicitly by 


py«'w = 2-“y 


+ a\ f n + 


k=0 


— k, 
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(x — l)^(x + 1) 


n—k 


Nevertheless, this formula is not the most useful one, at least for what we are looking 
for. Another equivalent characterization is the Rodrigues formula 


- 1 )-> + 1 )-'’ (I 


{x - (2.1) 


(cf. Section 4.3 of [63]). In particular, these expressions show that Jacobi polynomials 
are analytic functions of the parameters a,(3 G C and that degP^"’^^ < n. 

A third characterization of these polynomials, which will play the leading role in the 
sequel, is that they are the only polynomial solutions (up to a constant factor) of the 
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linear differential equation 


y"{x) + 


q; + 1^/3 + 1\ 

X — 1 X + 1 J 


y'{x) 



y{x) = 0, 


( 2 . 2 ) 


where A„ = n{n + a + j3 -\- 1). This is a second order linear differential equation of 
hypergeometric type. 


The setting is considered “classical” when a > —1 and (3 > —1. In this case, as we see 
from (2.2), the residues of the rational coefficient of y' in the differential equation are 
both positive. Moreover, we know that with this assumption the Jacobi polynomials 
are orthogonal on [—1,1] with respect the weight function (1 — x)°‘{l + x)^, i. e. 

J^^P^'^’f^\x)x’^{l-x)^{l + xfdx = 0, k = 0,l,...,n-l. (2.3) 

If both a > — 1 and (3 > —1, then the weight is integrable on [—1,1], and a well known 
consequence of (2.3) is that all the zeros of are simple and he on the interval 

(—1,1). Apparently, it was Stieltjes [61] who observed hrst that it is possible to give 
a nice interpretation of the location of these zeros as follows. Put two positive hxed 
charges of mass (/5 + l )/2 and ( q ; + 1)/2 at —1 and + 1 , respectively, and allow n positive 
unit charges X = {xi,... ,Xm} to move freely in (—1,1). If the interaction obeys the 
logarithmic potential law (that is, the force is inversely proportional to the relative 
distance), then in order to hnd the total energy E{X) of this system we have to add 
to the energy of the mutual interaction of these charges, 

A'mutual(A^) ^ ] 111 \Xk Xj \ , 

l</c<_7<n 

the component given by the “external held” (f{x) created by the hxed charges, 

(p{x) = |a; + 1| — —In |a; — 1|. (2.4) 

In other words, the total energy is 

n 

E{X) -^mutual (A) + E (p{xk) . (2.5) 

k=l 

There is a unique conhguration X* = {x\,... —1 < x]; < < ■ ■ ■ < x* < 1, 

providing the (strict) global minimum of E{X) in [—1,1]”, corresponding to the unique 
equilibrium position for our free charges. 


There uniqueness of the global minimum is not obvious, but there is an elegant proof 
in [63], based on the inequality between the arithmetic and geometric means. However, 
it is also a consequence of the following statement (Stieltjes’ theorem): points x\ are 
precisely the zeros of the polynomial The proof is quite straightforward and 

uses the diherential equation (2.2). Indeed, since E{X) +cx) as A approaches any 
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boundary point of [—1,1]"", we conclude that X* C (—1,1)"'. Then X* is also a critical 
point for the energy functional E{X), and as a consequence, every xl is a critical point 
of i? as a function of Xk (hxing the rest of a;*’s). Thus, the following necessary conditions 
must be satished: 


d 

dxk 


E{X) 


x=x* 


A. 

dxi 


E, 


mutual 


(x) 


x=x> 


+ ip'{xk) =0 for k 


1,2,...,n. (2.6) 


Let us consider the monic polynomial vanishing at x^’s: y{x) = Y\!j=i{.x — x*)-^ it is easy 
to check that 


d 

dxk 
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SO that (2.6) implies 


y"{x) — 2(p'{x)y'{x) = 0 for a; G X* . 


(2.7) 


For the external field (2.4) it means that 

y”{x) + (y'{x) = 0 for a; e X*. (2.8) 

\a; + 1 a; — 1 / 

However, the left hand side in (2.8) is a rational function of the type [n/2] with poles 
at ±1, so, up to a constant factor, it is equal to y{x)/{x^ — 1), which yields (2.2) and 
shows that y{x) = const pX^\x). 


Observe that we started from the strongest requirement of the global minimum of 
E{X) (that we refer to as a stable equilibrium) which implied that X* is also a point of 
local minimum, and in consequence, we have a “Nash-type” equilibrium (each function 
E{xk) = E(xl,..., xl_i, Xk, xl_^_i,..., x))^) attains its minimum at Xk = xl), yielding 
that each xl is a critical point of E{x\,..., xl_i, Xk, xIj^^, ... ,xD. Precisely this last 
statement turned out to be equivalent to the characterizing differential equation (2.2). 
The reader should take note of this clear hierarchy of equilibria in this problem. 


Interesting enough, in the situation we are considering all the equilibria are equivalent. 
In order to prove it we can consider the Hessian matrix 


H {hij) ) hij 


d^E{X) 

dxidxj 


of E{X). Straightforward computations show that 


(2.9) 


hij 


~2{xi - Xj) ^ if i 7^ j, 

+ 1 , CK + 1 1 ■ 

(a;* + l)2 {xi-1)^ ^ 


It is easy to observe that matrix E[ is real, symmetric, strictly diagonally dominant, and 
its diagonal entries are positive, from which we conclude that H is positive dehnite. In 
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particular, every critical point of E must be a point of a local minimum. Its uniqueness 
was observed above, hence it is in fact the point of the global minimum. 

Stieltjes considered also similar electrostatic models for other two classical orthogonal 
polynomials: Laguerre and Hermite. Since in this situation the free charges can move 
on an unbounded set, what can’t prevent them from escaping to inhnity? Stieltjes 
found a clever solution in putting a constraint on either the hrst (Laguerre) or second 
(Hermite) moment of their zero counting measures (see e.g. [66]). 

Namely, if for a > —1 we £x one charge of mass {a + l)/2 at the origin, and allow n 
positive unit charges X to move in [0, +cx)) with an additional constraint that their 
arithmetic mean is uniformly bounded, Y^^=iXi/n < K for some positive constant K, 
then again the unique configuration providing the global minimum of the total energy 
of the system of interacting particles coincides with the rescaled zeros of the Laguerre 
polynomial L^^^rnx), where r„ = (n + a)/K. In a similar way, allowing X to move 
freely on the whole R, but limiting the arithmetic mean of their squares, xj/n < K 

for some constant K, the unique configuration providing the global minimum of the 
total energy of the system of interacting particles coincides with the rescaled zeros of 
the Hermite polynomial Hn{snx), where = 

The proof for these both cases is similar to the Jacobi case, except that now we are 
dealing with a constrained minimum, so that the Lagrange multipliers corresponding 
to the restrictions when looking for a necessary condition of a critical point will become 
part of the characterizing differential equation. 

It is strange enough that Stieltjes himself did not realize that the unbounded cases 
(Laguerre and Hermite) could have been put in the same framework that the Jacobi 
one if we allow for an external held not necessarily generated by positive fixed charges. 
It was probably Ismail [30] who observed first that the zeros of these polynomials still 
provide the global minimum of the total energy (2.5) without any additional constraint 
or rescaling, if assuming for the Laguerre polynomials that (p is generated by a fixed 
charge of mass (a + l)/2 at the origin plus another term, equal to x/2; analogously, 
for the Hermite polynomials it is sufficient to take (p{x) = a;^/2. 

Reviewing the electrostatic models above several natural questions arise, such as: 

• are there generalizations of these models to other families of polynomials? 

• why necessarily the global minimum of the energy should be considered? Which other 
types of equilibria described above could be linked to the zeros of the polynomials? 

• what is the appropriate model for the complex zeros (when they exist)? 

• what kind of applications can these models have beyond their clear aesthetical value? 

In the following sections we will try to give at least some partial answers or to formulate 
conjectures related to these questions. 
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3 A generalization of the electrostatic model: Lame equation 


In fact, Stieltjes studied a more general situation. The generalized Lame differential 
equation (in an algebraic form) is 


E"ix) + E 


pi 


\i=0 


X 


E\x) 


C{x) 

A{x) 


E{x) = 0, A{x) = ^{x - ai), 


(3.1) 


i=0 


where C is a polynomial of degree < p — 1 (in the sequel we use the notation C G Pp-i). 
The case p = 1 corresponds to the hypergeometric differential equation, such as (2.2), 
whyle for p = 2 we obtain the Heun’s equation (see [51]). 

Heine [28] proved that for every n G N there exist at most 


a(n) 


fn + p-l\ 
[ n ) 


(3.2) 


different polynomials C in (3.1) such that this equation has a polynomial solution of 
degree n. These coefficients C are called Van Vleck polynomials and the corresponding 
polynomial solutions E are known as Heine-Stieltjes polynomials. 

Stieltjes studied the problem under the following two assumptions, generalizing the 
classical situation for the Jacobi polynomials: (i) the zeros of A are assumed to be 
simple and real, so without loss of generality we can take 


—1 = oo < oi < ■ ■ ■ < Op = 1, (3.3) 

(a) the residues of the coefficients of E' are all positive: 

Pi > 0, i = 0, ...,p. (3.4) 

(cf. equation (2.2) with assumptions a > —1, f3 > —1). Observe that zeros afs dehne 
p intervals {ok-i, ak) of the real line, and the combinatorial number a(n) dehned in 
(3.2) coincides with the number of different ways of placing n balls in p boxes. This 
is not a mere coincidence. According to Stieltjes [61], for each vector n = (ni,... ,np) 
such that ni + ■ ■ ■ + Up = n there exists a unique Heine-Stieltjes polynomial En and the 
unique corresponding Van Vleck polynomial Cn G Pp_i such that precisely zeros of 
En belong to the interval (afc_i,afc), k = 1,... ,p. Moreover, these zeros provide the 
unique global minimum for the total energy (2.5) under the additional restrietion of 
the number of zeros in each interval mentioned above, with the external held created 
by p + 1 positive charges hxed at afs: 


ip{x) 



(3.5) 
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This electrostatic model has been used in [41] in order to study the semiclassical limit 
n —> cx) (see Section 7), and in a number of papers [11,12] in describing the probabilistic 
distribution of these zeros as the number of intervals (that is, p) grows large together 
with n. 

Further generalizations of the work by Heine and Stieltjes followed several paths; we 
will mention only some of them. First, under assumptions (3.3)-(3.4) Van Vleck [68] 
and Bocher [9] proved that the zeros of C belong to [oq, Op]. A refinement of this result 
is due to several works of Shah [55,56,57,58]. Furthermore, Polya [50] showed that, 
allowing complex zeros of A, condition (3.4) is sufficient to assure that the zeros of 
E are located in the convex hull of a^’s. Marden [36] and later on, Alam, Al-Hashed, 
and Zaheer (see [1], [2], [69], [70]) established further results on location of the zeros 
of the Heine-Stieltjes polynomials under weaker conditions on the coefficients A and pi 
of (3.1). 

Going back to conhguration (3.3), not much can be said if we drop the condition of 
positivity of the residues pi in (3.4); the situation then becomes much more difficult to 
handle. The first attempt was made by Griinbaum [26,27], and later on by Dimitrov 
and Van Assche [18]. In particular, in [18] the case of p = 3 is analyzed when the 
positive and negative residues do not interlace: basically, the two cases considered are 
depicted in Figure 1. 


+ A+ - + A + 

• - • - • • •-• % -• 

Sq 3 i ^2 ^3 ^ ^1 ^2 ^3 


Fig. 1. Configuration studied by Dimitrov and Van Assche. 

Dimitrov and Van Assche proved that under certain restrictions on the degree of the 
polynomial and the residues pi, for every n E N there exists a unique pair (G„, En) of 
respectively Van Vleck and Heine-Stieltjes polynomials with degF^„ = n. All zeros of 
En belong to the interval enclosed by a/s with pj > 0 (denoted by A in Figure 1), 
and they are in the equilibrium position, given by the global minimum of the discrete 
energy (2.5) with the external held (3.5). As far as we know, no further studies in this 
direction exist. 

Interestingly enough, a similar situation appears when we study a special class of poly¬ 
nomials of multiple orthogonality (when the orthogonality conditions are distributed 
among several measures), which are object of an intensive study in the last years. 
Kaliaguine [32] (see also the more recent papers by Aptekarev [3] and Kaliaguine and 
Ronveaux [33]) considered the system of polynomials {p 2 n}, degp 2 n = 2n, satisfying 

J P2n{x)x^w{x) dx = J P2n{x)x’^w{x) dx = 0 , fc = 0, 1, . . . , U — 1 , 

with w{x) = x'^{l — x)"(l -|- x)^, ayjSyj > —1. It is a real case: p 2 n has exactly n 
zeros in each interval (—1,0) and (0,1). The asymptotic analysis (using the tools of 
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the potential theory) suggests that a reasonable model for these zeros might be that 
each of these sets of n zeros is the equilibrium position of the positive unit charges, 
which besides the mutual interaction are subject to the external field generated by the 
three hxed charges a/2, P/2, 7/2 at 1, —1, 0, respectively, and by the n charges of the 
weight 1/2 on the other subinterval. Nevertheless, this is an open problem. 

In a recent paper Grinshpan [25] considered a conhguration of zeros on the unit circle. 
He characterized the stable equilibrium of the n positive unit charges on T = G C : 
\z\ = 1} in the held generated by n negative unit charges in C\(Tn{0}) and showed that 
they are zeros of a Heine-Stieltjes polynomial satisfying a generalized Lame differential 
equation with coefficients symmetric with respect to T. Although it might give an 
impression that here we are dealing with complex zeros, the symmetry of the problem 
makes it essentially “real”. A more recent paper [40] extends the problem of stable 
equilibrium studied in [25] to the case of a held generated hj m = m{n) G N negative 
charges in C \ (T fl {0}) with values —uinki where oonk > 0 are the residues of Rn/Qn in 
the zeros Znk of Qn that he in the open unit disc, and Qn, Rn are polynomials of degree 
2m, 2m — 1 , respectively, such that, An{z) = zQn(z) and i?„ = Quiz) — zRn{z) are the 
coefficients of the corresponding generalized Lame diherential equation Aniz)y”iz) + 
Bn{z)y'{z) + Cn{z)y{z) = 0. 


4 Electrostatics and orthogonality 


Stieltjes’ theorem is beautiful, and also useful for proving some monotonicity properties 
of the zeros with respect to the parameters of the system. But is the existence of 
such a description a mere accident? Why should an electrostatic interpretation of the 
zeros of general orthogonal polynomials exist at all? One fact (known from the ’80) 
that makes us expect further developments is that the appropriate description of the 
asymptotic behavior (as n —> cx)) of the zeros of orthogonal polynomials is given in 
terms of some equilibrium measures, that is, measures minimizing certain logarithmic 
potentials, eventually in presence of an external held (see [53]). 

If we go back to the Ismail’s electrostatic model for the Laguerre and Hermite polyno¬ 
mials, then it is easy to recognize in the corresponding external held the inhuence of the 
weight of orthogonality of these polynomials. It was dehnitely an interesting problem to 
hnd out if this is a regular fact, tackled by Ismail in a series of papers [30,31]. We have 
noticed that so far the characterizing diherential equation has been the cornerstone of 
the model, and it was a reasonable line of attack. Thus, generalizing some pioneering 
works of Bauldry [7], and Bonan and Clark [10], Ismail deduced a second order linear 
diherential equation satished by a sequence of general orthogonal polynomials, after 
which he suggested a model for the zeros. Let us outline briehy his results. 

Let {Pn{x)} be polynomials orthonormal with respect to a unit weight function w{x) = 



exp(—v(a;)) supported on an interval [c,d] C M, finite or infinite: 


J Pm{x)pn{x)w{x) dx = 5m,n ■ 

Then sequence {pn\ satisfies a three-term recurrence relation 

XPn{x) = an+lPn+l{x) + hnPn{x) anPn-l{x), U > 0, > 0, 


with p-i{x) = 0 and Pq{x) = 1. Moreover, it has been shown in [7], [10], and also [17], 
[29], that under certain assumptions on w the orthonormal polynomials pn satisfy also 
the difference-differential relation 


p'^{x) = An{x)pn-l{x) - Bn{x)pn{x) , 

where the coefficients An and Bn are explicitly given in terms of tc, the recurrence 
coefficients a^, and the values of pn at the endpoints c and d. A direct consequence of 
this is that pn satishes also the second order linear differential equation 

Pn{x) - 2Rn{x)Pn{x) + Sn{x)pn{x) = 0 , (4.1) 


with 

^njx) 

2 ^ 2An{x) ’ 

Bn{x) - Bn{x) ^'^^^\ - Bn{x)[v'{x) + Bn{x)] + An{x)An-l{x) . 

-^n\X) dn—l 

In particular, if X* C is the ordered set of zeros of Pn, then 

Pnix) — 2Rn{x)p'n{x) =0 ioT X ^ X* . (4.2) 

Comparing it with (2.7) can associate (p'{x) = Rn{X), and it is natural to consider 
an electrostatic model for X* with logarithmic interaction between particles and an 
external held 

v{x) \n{knAn{x)) 

(p{x) = H- - -= y?iong(a:) + 9^short(a:) (4.3) 

{kn is any appropriate normalization constant, taken in [30] equal to Observe 

that this external held has two components: the hrst term in the right hand side of 
(4.3) has its origin in the orthogonality weight w{x) = exp(—n(a;)), and Ismail called 
it the long range potential. The second term, which received the name of the short 
range potential, is a bit mysterious, but explains several features of the classical models 
considered in Section 2, and at the same time, allows to give a further generalization 
of the electrostatic interpretation. 

The main result of Ismail [30] states that, assuming w{x) > 0 on (c, d), both v and 
ln(yl„) in C^(c, d), and the external held (4.3) convex, the total energy (2.5) has a 


Rn (^) 
Sn{x) = 
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unique point of global minimum, which is precisely X* (that is, the zeros of the or¬ 
thogonal polynomial Pn)- Obviously, what is used here is in fact the relation (4.2) 
equivalent to the fact that X* is just a critical point for E{X), and convexity assures 
that this is the unique point of minimum. The proof is in general very similar to the 
arguments for the Jacobi polynomials from Section 2. 

Beside some restrictions on the weight of orthogonality, one drawback of this model 
in application to general weights is that usually the short range potential cannot be 
explicitly computed in terms of w. This is not the case of the classical polynomials and 
some of their generalizations. 

For instance, for Jacobi polynomials, 

, , (1 — a;)“(l-f-a;)^r(Q;-|-+ 2 ) An{x) 2n + a + (3-\-l 

~ 2“+J+ir(a + l)r(/ 3 + 1) ’ an ~ 1 - ' 

In other words, w is responsible for the helds generated by two positive charges of 
mass {3/2 and a/2 at —1 and J-l, respectively, while the short range potential adds the 
missing charges of size 1/2 at ± 1 . 

Analogously, for Laguerre polynomials we have 

a;"exp(-a;) An{x) 1 

= -fw——^ , -= - , 

f (a -|- 1 ) an X 

so that the x°‘ factor in w plus the short range potential are responsible for the hxed 
charge of size (a + 1)/2 a.t z = 0 , and the exponential factor of the weight w generates 
the remaining component a;/2 of the external held ip, as described at the end of Section 

2 . 

So far the only role played by the short range component was adding some extra weight 
to the charges hxed by the weight. However, An can be responsible for a creation of 
“ghost” movable charges (either positive or negative) in the picture, as it was shown 
hrst by Grunbaum [26] (work that in fact predates [30]) or by Ismail himself when 
analyzing the case of the so-called Freud weight, w{x) = exp(—x^). 

Let us consider here a little bit more general situation. A generalized weight function 
(or measure) is semiclassical (see e.g. [37,59]) if it satishes the Pearson equation 

D{(f)w) = ipw , 

where 0, ip are polynomials, with degree of '^ > 1, and D is the “derivative” operator 
(in the usual, but also possibly in a distributional sense). It is well known that for such 
a weight the corresponding orthogonal polynomials (called also semielassieal) satisfy 
a diherential equation of the type (4.1), where the coefficients Rn and Sn are rational 
functions. The classical-type orthogonal polynomials considered by Grunbaum [26] are 
an example of a semiclassical family, but there are many more (for instance, the so- 
called sieved orthogonal polynomials [29], among others). 
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The systematic study of the semiclassical orthogonal polynomials from the point of view 
of the electrostatic interpretation of their zeros has started in the works of Garrido and 
collaborators [23], [24], For instance, they consider a perturbation of the Freud weight 
function {w{x) = exp(—by the addition of a fixed charged point of mass A at the 
origin; the corresponding orthogonal polynomials are known as Freud-type polynomials. 
The resulting orthogonality measure is semiclassical, and it was proved in [23] that these 
polynomials satisfy a second order linear differential equation of the form (4.1), and 
the electrostatic model is in sight. In fact, in the situation considered, and following 
Ismail’s terminology, the long range potential is, as expected, 9?iong(a^) = a;^/2, while 
the short range potential depends on the degree n and has the following structure: 


9-^ short ( 2 ^) 


1 Qi{x,n) 

2 


(4.4) 


where Qi{x,n) has two real roots ri(n),r 2 (n) and two simple conjugate complex roots 
r 3 (n),r 4 (n). Their asymptotic behavior is 

hmri(n) = limr 2 (n) = 0 , hmr 3 (n) = limr 4 (n) = 00 . (4.5) 

In fact, in the odd case, both real zeros coalesce at the origin: 


ri(2n -I- 1) = r2(2n -|- 1) = 0 , 

canceling out the denominator in the argument of the logarithm in (4.4). In other words, 
the short range potential is given by one unit positive charge fixed at the origin, plus 4 
floating negative charges of weight —1/2 on C situated symmetrically with respect to 
the origin. However, in the odd case two of these negative charges annihilate with the 
fixed positive one, paving the way to a zero of the orthogonal polynomial to take the 
origin, position it must occupy due to the symmetry of the problem. This annihilation 
has as a consequence that for odd n’s the Freud and the Freud-type polynomials are 
identical. Moreover, (4.5) shows that also for the even case they are asymptotically 
identical (as n cx)). 

One of the main results of [23], [24] is that the zeros of the Freud-type polynomials 
provide a critical configuration for the total energy E{X) (in other words, that for 
these zeros X* equations (2.6) are satished). But can we assure that they are in a 
stable equilibrium, that is, that E{X*) = mini?(X)? The complete answer is not clear 
yet. By compnting the Hessian matrix H (2.9) it was proved that for small values of 
the correction charge A this is really the case, bnt there are some evidences that make 
it possible to conjecture the existence of a critical valne Aq snch that for 0 < A < Ao(n) 
the zeros of the orthogonal polynomials attain the minimum of the total energy, but 
for A > Ao(n) the equilibrium is no longer stable. Is still any other type of equilibria 
discussed in Section 2 preserved in this case (beside the critical point of the energy)? 
This is also an open question. 
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5 A mcix-min problem 


From potential theory it is well known that many of energy minimization problems 
have in fact a min-max (or max-min) characterization. Let us analyze again, but from 
this point of view, the classical models of Stieltjes, and in passing, address the issue 
why the real line should be present in the model. 

This question is not trivial; let us go back to the Jacobi polynomials with a, /3 > 

—1. Why their zeros should be real? A standard explanation is that these polynomials 
satisfy the orthogonality conditions (2.3), but where the real line is present in that 
integral? Observe that the integrand in (2.3) is an analytic function, so we perfectly 
can obtain the same result integrating along any reasonable curve joining —1 with + 1 . 

The real line is also built in the electrostatic model of Stieltjes as an a priori constraint. 
However, we may put forward an alternative model free from this restriction. As in the 
Stieltjes’ setting we fix two positive charges of mass (/? + l )/2 and (a + l)/2 at —1 
and +1, respectively, and denote again by E{X) the total energy of a discrete system 
X C C”, given by formulas (2.4)-(2.5) (assuming that in case of coincidence of any two 
points, E{X) = +CX)). Denote by T the family of compact continua on C containing 
both —1 and +1. For any A G JF set 

m{K) = M{E{X) : X = {x^ ■ ■ ■ ,xn} C K} , 

and consider the following extremal problem: 

E* = sup{m(A) : K E E} . 


Theorem 1 E* < +cx), and this value is attained at a unique n-tuple X* = {xj;, ..., x*}, 
characterized by the fact that x* are precisely the zeros of the Jacobi polynomial 

Proof. Let X be any conhguration of n distinct points on the interval (—1,1) (the 
case of coalescence of some points is trivial), and let A C JF be any other continuum. 
Since ±1 G A, we can assure that there exists a discrete set of points Z C A of the 
form 

Z = {zj = Xj + iyj : G M, j = 1,..., n} 

(that is, any vertical line passing through an Xj must intersect A). It is straightforward 
to check that 


\xk — Xj\ < \zk — Zj\ , k ^ j , and |1 ± < |1 ± 2:^1 , j = 1,. .., n. 

Taking into account the expression of the total energy and the positivity of all the 
charges involved in the system we immediately conclude that E{X) > E{Z). Hence, 

m([—1,1]) > m{K ), 
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and since K was arbitrary, we obtain that E* = m([—1,1]). But from Stieltjes’ theorem 
it follows that there exists a unique X* C (—1,1)” such that E{X*) = m([—1,1]), and 
the points of X* are the zeros of P^’^\ which concludes the proof. □ 

Obviously, a similar result is possible to formulate for Laguerre and Hermite polyno¬ 
mials. 

The value of the min-max model we have described is that it allows to consider easily 
complex zeros of these classical families, which appear when the parameters involved 
become no longer “classical” (see for instance Figure 2 with two examples of zeros 
of Jacobi polynomials). In fact, as it follows from the Rodrigues formula (2.1), if we 
allow for arbitrary real a and (3, the Jacobi polynomials constitute a very rich object 
exhibiting several types of orthogonalities: hermitian, non-hermitian, and even multiple, 
depending on the values of the parameters and the degree (see the classihcation of all 
cases obtained [34]). We conjecture that Theorem 1 is valid in all these situations after 
a suitable modification of the family of continua T. For instance, ii 2n + a + j3 < —1 
(see Figure 2, left), the appropriate class T would include all continua extending to 
inhnity, and separating —1 and J-l. From the results of [34] a rule of the thumb follows: 
if a + n > 0 (resp., if (3 + n > 0) then the continua from T must contain J-l (resp., —1). 
Moreover, ifn-|-Q;-|-/3-|-l<0, then the continua from T must include the inhnity. 



All these facts have been conjectured by one of this authors some time ago, but their 
proof still remains an open problem. A general approach to such kind of theorems 
would allow also to tackle other kind of problems not linked directly to the real line or 
exhibiting non-real zeros. 

Typical examples are those related to the multiple (or Hermite-Pade) orthogonality, 
when the orthogonality conditions are distributed among several measures, or of non- 
hermitian orthogonality, when the path of integration is, in fact, not hxed a priori. In 
many situations these polynomials exhibit complex zeros, whose asymptotic distribu¬ 
tion has been studied by several authors. From the pioneering works of Stahl [60] it is 
known that the limiting location of the zeros is usually described by some trajectories 
of a quadratic differential, which in turn have a min-max description. This gives an ad¬ 
ditional evidence to our conjecture, that in fact says that this behavior is not exclusive 
of large n’s (a kind of “discrete version” of Stahl’s theorem). 
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6 Further generalizations and applications of the electrostatic model 


Many new interesting problems related to the electrostatic models have their origin 
either in some non-standard orthogonality conditions (such as non-hermitian or multi¬ 
ple orthogonality mentioned above), or in physical applications. In this Section we will 
briefly review some of them. 

Exactly solvable or quasi-exactly solvable multi-particle quantum mechanical systems 
have many remarkable properties. Especially, those of the Calogero-Sutherland-Moser 
(CSM) systems [13,14,45,62] and their integrable deformation called the Ruijsenaars- 
Schneider-van Diejen (RSvD) systems [52,67] have been well studied. A classical result 
is that the equilibrium positions of the CSM systems are described by the zeros of 
the classical orthogonal polynomials; the Hermite, Laguerre, and Jacobi polynomials 
[15,16] (see also [48] for a comprehensive review and bibliography therein). Following 
this analogy the authors of [48] proved recently that the equilibrium positions of the 
RSvD systems with rational and trigonometric potentials coincide again with zeros of 
polynomials from the Askey tableau of hypergeometric orthogonal polynomials; namely, 
they found connections with the Meixner-Pollaczek, continuous Hahn, continuous dual 
Hahn, and the Askey-Wilson polynomials. Hence, an electrostatic model for the zeros 
of hypergeometric polynomials in all their generality could have several interesting 
applications for exactly solvable quantum-mechanical systems. 

Another interesting model was considered by Loutsenko [35], who studied a system 
(X, Y) C of n positive and m negative moving charges in C of masses 1 and 

—A, respectively (with A G {1/2,1, 2}), interacting again according to the logarithmic 
law. He proved in particular that if {X*,Y*) is a critical configuration of the mutual 
energy i7jnutuai(Ai, F), p is the monic polynomial vanishing at the points of X*, and q 
is the monic polynomial vanishing at those of Y*, then p,q are solutions of a bilinear 
differential equation 

p"q — 2kp'q' + k^pq” = 0. 

Regarding the Lame equation, once we depart from the classical setting considered by 
Stieltjes, the description of the Heine-Stieltjes polynomials is not apparent. It is clear 
however that we can hardly expect that their zeros provide a stable equilibrium; instead, 
we must concentrate on the least restrictive condition of configuration providing critical 
values of the energy functional. This study has again some physical applications. The 
BCn elliptic Inozemtsev model is a quantum integrable system with n-particles whose 
potential is given by elliptic functions. For the case n = 1, finding eigenstates of its 
Hamiltonian is equivalent to solving the Heun equation (see e.g. [64,65]. In this sense, 
the BCn Inozemtsev model is a generalization of the Heun equation. 

There are also connections of some problems in physics and representation theory with 
multiple orthogonal polynomials. It is known that zeros of the Jacobi polynomial 
satisfy a system of algebraic equations, which is known as the Bethe Ansatz equation of 
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the Gaudin model associated to s /2 and two irreducible modules with highest weights 
— (a+l), —(/3+1). In [47] the authors generalized this connection, studying sequences of 
r polynomials whose zeros constitute the unique solution of the Bethe Ansatz equation 
associated with two highest weight slr+i irreducible modules, with the restriction that 
the highest weight of one of the modules is a multiple of the first fundamental weight. 
As a result, they show that the hrst polynomial in the sequence coincides with the 
well known Jacobi-Piheiro multiple orthogonal polynomial, and others are given by 
Wronskian type determinants of Jacobi-Piheiro polynomials. In [46] they derived a 
linear differential equation (alas, of order r) for these polynomials, paving a way to the 
electrostatic interpretation of their zeros. 

Finally, we can mention the matrix orthogonal polynomials, whose non-trivial connec¬ 
tions with (matrix) differential equations is being studied. Duran and Griinbaum gave 
in [19,20] some examples of matrix orthogonal polynomials {Pn} satisfying the second 
order linear ODE of the form 

P”{x)A 2 {x) + P^{x)Ai{x) + Pnix)Ao{x) = TnPnix), U > 0 , 

where (as in the scalar hypergeometric case) coefficients Aj are matrix polynomials 
that do not depend on n, of degrees < 2, 1, and 0, respectively, and P^ are Hermitian 
matrices. These {Pn} look like a natural generalization of the families of classical poly¬ 
nomials in the scalar case; likewise, they exhibit a rich variety of structural properties 
(see [21]). Furthermore, they might have some relevance in the analysis of the Dirac 
equation (relativistic analogue of the Schrodinger equation). However, a study of the 
electrostatic interpretation of the zeros of the orthogonal matrix polynomials remains 
completely open. 


7 Asymptotic distribution of zeros 


Assume we have an electrostatic model of n positive charges moving in an external 
field, and we are interested to analyze what happens when n ^ oo (the so-called 
thermodynamic or semiclassical limit). In many situations this information can be 
extracted directly from the characterizing second order linear differential equation, 
using either the WKB approach (cf. [4]-[6], [42]-[44], [71]-[74]), or reducing the ODE 
to a Riccati form, a method that we outline below. But first we should introduce some 
notation and agree in the meaning of the “global” asymptotics of the zeros. 

In the sequel, supp(/i) denotes the support of a measure /i, and 

^ f z eC\ supp(/r), 

J z — t 

is its Stieltjes (or Gauchy) transform. With a function y : C ^ C we can associate its 
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zero counting measure Uy, defined by 


y{x)=0 


where S^; is the unit mass (Dirac delta) at x, and the sum goes along all the zeros of 
y taking into account their multiplicity. Equivalently, for each Borel set A C C the 
number of zeros of y in A is 

z/j,(A) = / duy{x) . 

JA 

An important observation is that if p is a polynomial and Up is the associated zero¬ 
counting measure, then the logarithmic derivative of p is the Stieltjes transform of v. 


P\z) ^ 1 

p{z) ^^z-ak 


du{t) 
z — t 


= 1Z[Z 


which can be evaluated away from the zeros of p. 


(7.1) 


The “global” behavior of the zeros of polynomials {p„}, deg(p„) = n, is described by 
the limit of the sequence of normalized measures pn = iZp^/n as n ^ oo, in the sense of 
the weak-* topology. Recall (cf. [8, Section 2]) that a sequence of Borel measures {pn} 
converges to a measure p in the weak-* topology (which we denote as pn —^ p) if 

, V/ continuous and compactly supported on C , 

Another important fact is that the set of unit measures with uniformly bounded sup¬ 
ports is compact in the weak-* topology. Hence, if we know for instance that all the 
zeros of the sequence of polynomials {pn}, deg(p„) = n, belong to the same compact 
set K C C, then there always exists a unit measure p supported on K and a subse¬ 
quence A C N such that /i„ ^ p for n G A (where pn are the normalized zero counting 
measures dehned above). In consequence, and taking into account (7.1), we have that 


lim — = jgl^z ), 2 : G C \ supp(p). (7.2) 

nGA n Pn[z) 

The method of reduction of the original second order ODE to the Riccati form is based 
on the previous observation. This idea appears in a work of Saff, Ullman, and Varga 
[54], although its roots can be traced back to the famous Perron’s monograph [49], who 
in turn gives credit to some original works of Euler. Recently it has been successfully 
applied in a variety of problems (see, e. g. [22,38,39,44]). 

Let us see how it works in the simplest case of Jacobi polynomials p„ = with 

a,P > —1. The first step is to rewrite equation (2.2) in terms of the normalized 
logarithmic derivative of the polynomial solution. 


hn{z) 


1 Pnjz) 
n Pn{z) ’ 
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which yields 


'K{^) 


n 


+ h\{z) + 


[x\ 


Cl + l^/3 + l\ hn ^ 
z — 1 z + 1 ) n 


+ n + a + )3+1 ^0. (7.3) 


n 


Since all the zeros of pn are in [—1,1], by the argument explained above we can assure 
the existence of a unit measure p, supported on [—1,1] and a subsequence A C N such 
that (see (7.2)) hn{z) p{z), n G A, uniformly on compact subsets of C \ [—1,1], 
where h'^ are also uniformly bounded. Taking limits in (7.3) when n G A we arrive at 
an algebraic equation for h = p\ (1 — x^)h?'{x) + 1 = 0, so that 

aW = 


and the appropriate branch of the square root in C \ [—1,1] is fixed by the condition 
Wmz^oo zp{z) = 1. Furthermore, since p is independent on A, we can claim that p is 
in fact the weak-* limit of {pn} along n oo. 


It remains only to recover p from its Stieltjes transform. Since we know (important 
fact!) that supp(/i) C [—1,1], this task is just a straightforward application of the 
Sokhotsky-Plemelj’s formulas (that allow to find p from the boundary values of p on 
(—1,1)), that yield that the limit measure p is absolutely continuous and its density is 

/ / N 1 dX ^ , N 

= — /. 2 > 0 ^ a;G(-l,l), 

TT V1 — 

(fact that we do not claim that was discovered by us; probably, it is at least 150 years 
old). 


In some more difficult situations the described elementary method becomes technically 
involved, and we have to turn to other resources. One of the advantages of the elec¬ 
trostatic interpretation of the zeros of a sequence of polynomials is that it allows us 
to guess what happens in the thermodynamic limit. For instance, if for each n the 
configuration of zeros is a global minimizer of the total energy, then it is natural to 
expect (and can be proved) that the asymptotics in a “global sense” should be well 
described (at least, in the first approximation) by a continuous distribution which is 
solution of the corresponding extremal problem for the logarithmic potential energy in 
the class of all probability measures (including the absolutely continuous ones). This 
distribution is known as the equilibrium measure, probably, in an external field (see 
[53] for details and definitions). 

Let us consider for example the generalized Lame equation (3.1) under the assumptions 
(3.3)-(3.4). We are interested in the limit n ^ oo with the extra assumptions 

lim — = 6'i, i = 1 ,... ,p 

n—>oo ^ 

(see Section 3 for notation). Since for every n the zeros minimize the total discrete 
energy of the system, we can expect that any weak-* limit p of the zero-counting 
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unit measures of the Heine-Stieltjes polynomials En solves a similar, but continuous, 
extremal problem. Indeed, it was proved in [41] that \i M. = M(0i,..., dp) denotes the 
class of all probability measures i/ on [— 1 , 1 ] such that 



di, i = 


then /i is the minimizer of the logarithmic energy in the class J\4. Existence and unique¬ 
ness of this minimizer is proved by standard methods of potential theory. At this stage 
we can apply the reduction to the Riccati form method described above to get the fnll 
description of p. Let us summarize briefly some results from [41]. 


For any system of p — 1 points 


-1 < A < ■ ■ ■ < /9p-i < 1 


(7.4) 


we dehne the fnnctions 


R{x) := [Q (x - pj ), H{x) := ^ 


R{x) 

A{x) 


-, X 
X 


oo. 


Then there exist p — 1 points (7.4) nniqnely determined by the following system of 
equations: 

Im / H{x) dx = —TT dj , j = 1,..., p — 1, 

Jaj-i 

where we take the limit values of H from the upper half plane. If we introdnce the 
counting function 

Z{x) := [va - vr] ^(-oo,a;]^ , 

then snpp(/i) = {x G M : Z{x) = 1} and it consists of at most p — 1 disjoint intervals 
in [-1,1]. 


Furthermore, p is an absolutely continuous measure, 

1 1 

p'{x) = - H{x) = — \H{x)\ , X G supp(/i), 

m n 

and, for 2 : ^ snpp(p), p{z) = —H{z). 


The methods described in this section yield many nice results, but can fail for two 
reasons: 


(a) If we have no a priori information on the location of the zeros, then the reduc¬ 
tion of the ODE to the Riccati form will give us at most an expression of the 
Cauchy transform p of the limit distribution p in the domains (nnknown) disjoint 
with snpp(/i) (and supp(p) could be, eventually, a subset of C of positive plane 
measnre). Does this information determine pi It is not clear, althongh for some 
specihc expressions of p this should be really the case. 
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(b) If we are analyzing Heine-Stieltjes polynomials whose zeros provide a critical con¬ 
figuration for the total energy E{X) (but not the global minimum!), then what 
measure p should we expect in the semiclassical limit? Clearly, p is not necessarily 
an equilibrium measure, but rather a critical measure, that on the real line can 
be characterized by the fact that its potential (plus the external held, if exists) is 
constant on each connected component of supp(p), but unlike in the equilibrium 
case, these constants need not to be the same. However, in this description two 
open problems remain: (i) io prove that discrete critical measures converge in 
the thermodynamic limit to continuous critical measures, and that all continuous 
critical measures can be obtained this way; (ii) to hnd a feasible description of 
the multi-parametric family of continuous critical measures in the given class. 
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